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Abstract 

A Faraday rotation experiment can set limits on the magnetic moment of a electrically-neutral, 
dark-matter particle, and the limits increase in stringency as the candidate-particle mass decreases. 
Consequently, if we assume the dark-matter particle to be a thermal relic, our most stringent 
constraints emerge at the keV mass scale. We discuss how such an experiment could be realized and 
determine the limits on the magnetic moment as a function of mass which follow given demonstrated 
experimental capacities. 
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I. INTRODUCTION 



Disparate astronomical observations provide compelling evidence for additional, non- 
luminous matter, or dark matter, in gravitational interactions. The evidence includes the 
persistence of the galactic rotation curves to distances for which little luminous matter is 
present [H E], the relative strength and shape of the galaxy-distribution power spectrum at 
large wave numbers [3j, and the pattern of acoustic oscillations in the power spectrum of 
the cosmic microwave background [3]. The cosmological evidence, in aggregrate, points to 
the assessment that dark matter comprises some twenty-three percent of the energy density 
of the universe, with a precision of a few percent [5J. Yet we know little about the nature 
of dark matter — we do not know its mass, its quantum numbers, or even with surety that 
it is indeed composed of isolated elementary particles. A recent gravitational lensing study 
does disfavor a modification of gravity in explanation of its effects [S]: we shall assume [7| 
that dark matter exists and is matter. 

We do know other things about dark matter [HI [9]; namely, that it is not hot [10], and 
that it appears to lack both electric and color charge [H] [121 [El [El [iS, [16] . Here we use 
temperature, i.e., whether dark matter is "cold" or "hot," to connote whether dark matter is 
non-relativistic or relativistic, respectively, at the redshift at which it decouples from matter 
in the cooling early Universe. For so-called thermal relics, this criterion selects the mass 
of the dark candidate as well, so that colder particles are heavier. However, alternative 
production scenarios can exist, and very light particles can also act as cold dark matter, 
as in the case of the axion |18j. Nevertheless, model-dependent constraints do exist on the 
mass of the candidate particle, and we consider them in Sec. |Tl} 

The evidence for dark matter emerges from astrophysical observations of gravitational 
interactions, but establishing its couphngs to Standard Model (SM) particles has proven 
elusive. Nevertheless, such a quest is of great import, for it is through such means that its 
mass and quantum numbers can ultimately be determined. Indeed all direct and indirect 
means of detecting dark matter rely on the notion that it does indeed experience weak or 
electromagnetic interactions to some degree. Turning to the known particles of the SM for 
guidance, the neutron shows us that a particle can have both a vanishing electric charge and 
a significant magnetic moment. Thus we wish to constrain the possibility that dark matter 
has a small electromagnetic coupling, via its magnetic moment. We review the existing 
constraints on this possibility in Sec. |III[ 

In this paper we consider a new technique for the direct detection of dark matter, namely, 
through use of the gyromagnetic Faraday effect [19]. Alternatively, this effect can be used 
to limit a possible magnetic moment /i of a dark-matter particle of mass M. Let us consider 
how this could work. An electrically neutral medium of particles which possesses a net 
magnetization in an external magnetic field is circularly birefringent, even if the medium is 
isotropic. This implies that the propagation speed of light in the medium depends on the 
state of its circular polarization, so that light prepared in a state of linear polarization will 
suffer a rotation of the plane of that polarization upon transmission through the medium [20l 
[21] , as long as it does not travel at right angles to the external magnetic field. We term this 
the gyromagnetic Faraday effect, after Ref. |22]. Note that we need not rely on any existing 
magnetization of the dark matter in order to realize an effect. Rather we imagine a Faraday 
rotation experiment mounted in a region with a large external magnetic induction Bq and 
dark matter of spin 5* incident on it in the direction of Bq. If the value of iiBq is larger 
than the dark matter particle's kinetic energy in the Earth's rest frame, then the field region 
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acts as a spin-filter, or longitudinal Stern-Gerlach, device — at least the highest energy spin 
configuration cannot enter the field region. This technique is used to polarize ultra-cold 
neutrons (UCNs) in the UCNA experiment at Los Alamos with near 100% efficiency 
Thus, viewed in the Earth's rest frame, the dark matter which sweeps through the field 
region can possess a net magnetization. If light transits this medium in the direction of 
the magnetic field, and if we define k± to be the wave number for states with right- (+) or 
left-handed (— ) circular polarization, then the rotation linearly polarized light suffers in its 
transit through the medium is given by the angle (f) = — /c_)//2, where / is the length of 
transmission through the medium. If (p is non-zero once all systematic effects which could 
mimic the signal are excluded, we have evidence for dark matter with a non-zero magnetic 
moment. 

The direct detection of dark matter with a magnetic moment could be realized in a 
variety of ways. For example, one could search for anomalous recoil events in scattering 
from nuclei, in just the manner one searches for spin-independent and spin-dependent dark 
matter-nucleon interactions [211 iSSl [26] . Indeed, an experimental signal of the latter, from 
the DAMA experiment [27], has been interpreted as a limit on a putative dark- matter 
magnetic moment, namely of /x < 1.4 • 10~^/iAr for a dark matter candidate of 100 GeV in 
mass [28]. In such experiments, however, the light dark-matter candidates we discuss give 
rise to momentum transfers which are much too small to be detected, even if the channeling 
effect proposed by Ref. and studied in Ref. is operative for the events studied in 
the DAMA/LIBRA Nal detector [31J — note Ref. [32 for further discussion. Alternatively, 
a dark matter magnetic moment can be found through either nuclear magnetic resonance 
(NMR) or Faraday rotation studies. A NMR signal is typically realized through the detected 
change in magnetic field incurred through a spin-fiip transition of the magnetic-moment- 
carrying particle, induced by an applied radio frequency. However, very small magnetic 
fields are more efficiently discovered through magnetoptical studies [331 [Mj. Thus we focus 
on the use of Faraday rotation to detect dark matter. Precision optical rotation studies 
have also been conducted in which the external magnetic field is oriented at right angles to 
the direction of the propagating light, to the end, e.g., of testing the optical birefringence 
of the vacuum, or, alternatively, of limiting the photon-axion coupling ^35] EB]. Such an 
experimental configuration does not support a nonzero Faraday effect, yet the empirical 
parameters used in these studies are useful to us, for we employ them to estimate the limit 
on the magnetic moment as a function of mass in the set-up we consider. 

A great variety of dark matter candidates, consistent with the various astrophysical con- 
straints, exist in the literature [371 EH]- Indeed, their masses vary from some 10~^^ to 10^^ 
GeV, and their interaction cross sections — with nucleons — vary over orders of magnitude 
as well. Typical direct detection strategies rely on the observation of "anomalous" nuclear 
recoils [2S1 12S1 EH], so that their sensitivity is typically to candidates of C(100 GeV) in mass 
scale. The experiment we suggest is sensitive to a completely different window in param- 
eter space — to dark-matter candidates of crudely 0(1 MeV) or less in mass. Although 
constraints can be set through cosmological studies [I9] , terrestrial studies are amenable to 
better control, for the existence of cosmological magnetic fields have not yet been estab- 
lished [39]. 

Let us conclude our introduction by outlining the sections to follow. We begin by de- 
scribing, in Sees. [IT] and III, respectively, the mass and magnetic moment constraints which 
exist on a light dark-matter candidate with a non-zero magnetic moment. We then proceed 
to review the gyromagnetic Faraday effect in Sec. |IV| and to describe in concrete terms the 
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experimental limits on /i one might possibly attain in Sec. |Vj We conclude with a summary 
in Sec. ED 

II. MASS CONTRAINTS 

In standard Big-Bang cosmology, the nature of dark matter impacts the formation of the 
large-scale structure of the Universe. In particular, if dark matter is cold and coUisionless, 
then galaxy formation proceeds via a hierarchial clustering, namely, from the merging of 
small protogalactic clumps, on ever larger scales [101 HH 1121 US] • In contrast, if dark matter 
is hot, the hierarchy is inverted, so that large protogalactic disks form first, which then 
clump im US]. Galaxies, however, are observed at much larger redshifts than such simu- 
lations predict |H1 Il2]. Moreover, observations of particular classes of quasar absorption 
lines, the so-called damped Lyman-a systems, thought to be the evolutionary progenitors 
of galaxies today, also favor the former scenario pH]. It has also been argued that hot dark 
matter, i.e., most notably, light, massive neutrinos, cannot explain the galactic rotation 
curves [15]. However, the cold-dark-matter paradigm does have difficulties in confronting 
small-scale structure; it yields, in effect, too much dumpiness below the Mpc scale. Warm 
dark matter has been advocated as a way to alleviate these difficulties Limits on the 
mass of warm dark matter emerge from the comparison of the observations of the Lyman-a 
absorption spectrum with numerical simulations [l?! HHl HHl ED] ; the limits depend on the 
particle considered, and the manner in which it is produced [51], yielding [50], at 2(T, M > 4 
keV for a thermal relic and M > 28 keV for a massive sterile neutrino [5^ . 

Cosmological constraints also exist on the mass of a dark-matter particle. If the particles 
annihilate via the weak interaction, then aannV is parametrically set by JVaC^^^M^ , where 
Gp is the Fermi constant. Ma is a dimensionless factor, and we assume aann oc 1/v. In 
this case avoiding a dark-matter abundance in excess of the observed relic density bounds 
M from below. Indeed, under these conditions the mass of the cold dark-matter particle 
must exceed (9(2 GeV) to avoid closing the Universe [53]. The resulting lower bound on 
M can be relaxed in different ways. Feng and Kumar [5l], e.g., have emphasized that the 
appearance of Gp in (TannV is simply parametric, that Gp can be replaced with g^s, and that 
the effective coupling Qes can be small without having the precise numerical value of Gp- 
Thus if Qes > Gp, the bound on M is weakened. Indeed, such considerations permit dark 
matter candidates which confront the relic density and big-bang nucleosynthesis constraints 
successfully but yet range from the keV to the TeV scale in mass [HH [SS] . 

In this paper we consider what empirical constraints can be placed on the magnetic mo- 
ment of a dark matter particle. This hypothesis gives rise to a new annihilation mechanism, 
though both a dark matter particle and its antiparticle must be present to realize it. We 
recall that particles with magnetic moments are invariably described by complex field rep- 
resentations, so that such a particle and its anti-particle are physically distinct — by the 
CPT theorem we expect the magnetic moments of such particles to differ only in sign. If the 
particle-antiparticle annihilation is mediated by a magnetic moment interaction, then aannV 
is parametrically set by M'^a"^ fj,'^ with Af^ a dimensionless parameter, as long as M > me, 
the electron mass. The annihilation of still lighter mass dark matter candidates follows a 
different parametric form. That is, if the particles annihilate to Standard Model particles, 
then o"a„„f is suppressed by higher powers of the coupling constant — at least. Alternatively, 
if they annihilate to "secluded" dark matter particles, note, e.g., Ref. [56], then aannV is of 
form A/'jiaa'/x^, where a' is the electromagnetic coupling of the secluded particles and Af^ 
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is a dimensionless parameter. Generally, we expect the magnetic moment /i to be of form 
II = Keh/2M, where k is the anomalous magnetic moment, so that (TannV scales as 
The presence of an additional annihilation mechanism should make the light dark matter 
candidates we consider, of C(keV) scale in mass, say, decouple as matter and not radiation, 
so that constraints on the number of relativistic degrees of freedom during the epoch of big- 
bang nucleosynthesis also do not apply. This new annihilation mechanism becomes more 
effective as the candidate mass grows lighter. Nevertheless it is still possible to saturate 
the dark-matter density with such a candidate particle, for the efficacy of the annihilation 
process can be mitigated by a particle-antiparticle excess |571 |58]. Thus the dark matter 
relic density need not bound the magnetic moment from above. Moreover, we emphasize 
that dark matter could have multiple components, so that an upper bound on the magnetic 
moment could also be evaded by diluting the magnetic-moment-carrying particle with other 
sorts of dark matter. We shall assume these various annihilation mechanisms are effective 
enough to permit dark matter candidates as light as (9(1 keV) in mass. 

The intensity and morphology of galactic positron emission, as studied by the INTEGRAL 
satellite |59j , has prompted much discussion of dark matter candidates with electromagnetic 
interactions |5H |60l EB EH |63l EU EHl ES] , which are of 0(1 MeV) scale in mass, as well as of 
other possibilities, as, e.g., in Refs. [671 ES ESI [70] . If the pattern of the INTEGRAL spectra 
are indeed explained by dark matter, then additional constraints follow on its nature. For 
example [nil E2] , observational constraints on the diffuse photon flux also impose limits on 
the mass of the dark-matter candidate x through internal bremsstrahlung corrections to 
the annihilation process XX ~^ e~^e~. Such constraints apply to our scenario as well; the 
upshot is that the dark particle's mass is limited to be less than a few MeV [6T1 [62]. In our 
case internal bremsstrahlung contributions can also be generated by the magnetic-moment- 
carrying particle. However, we note that soft photon emission via a Ml transition of an 
isolated magnetic dipole is slow compared to the rate set by the inverse age of the universe, 
as we discuss in greater detail in Sec. [IV[ so that no meaningful limit follows on its magnetic 
moment. In what follows we consider candidate particles which range from 0{1 keV) to 
(9(100 MeV) in mass, though our limits are most effective at sub-MeV mass scales. 

III. MAGNETIC MOMENT CONTRAINTS 

Various constraints on the magnetic moment of a dark-matter particle for masses in 
excess of 1 MeV have been considered in Ref. [SB]. We review these and more, in order 
to provide a context for the direct detection experiment we suggest. In the mass window 
of interest to us, two experimental constraints are important — one comes from precision 
electroweak measurements [58], and the other comes from low-energy e"'"e~ collider data, 
namely from the process e~^e~ —* vv^ [ZH [72]. This last constrains the magnetic moment of 
the invisible particle directly and thus offers a constraint on a dark-matter magnetic moment 
as well. The authors of these studies use data at center-of-mass energies sufficient to produce 
a Ut-Ut pair, given accelerator constraints on its mass [71] — this easily includes our mass 
range of interest. Interpreting their results as a limit on the anomalous magnetic moment 
of the tau neutrino, the low-energy analyses conclude < 4 ■ 10~^/iB at 90% CL [7T] . 
and iiy^ < 9.1 • 10"^/^^ at 90% CL [72] from distinct data sets. A more severe limit on the 
Vt magnetic moment does exist [75]; however, the nature of the e^e~ limits allows us to 
interpret them in a manner useful to our current study. For a discussion of how low-energy 
e^e~ collider data can probe particular MeV dark matter models [60], see Ref. |74j . 
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Precision electroweak measurements also constrain the magnetic moment [58] • The quan- 
tity Af captures radiative corrections to the relationship between the fine-structure constant 
a, the Fermi constant Gp, and the and Z masses, Mw and Mz [73]. The differ- 
ence between the empirically determined value of Af and that computed in the Standard 
Model provides a window Af"^*^™ to which a dark-matter particle can contribute. Following 
Refs. [HHl [76] , we assume Af°°" is given by the vacuum polarization correction to the photon 
self-energy from a dark-matter particle with a magnetic moment, with no other adjustments. 
We choose to study the quantity Af as its uncertainty is dominated by that in the running 
of a |T7j . Thus we consider [75] 
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where Sz is computed in the MS scheme and is (1 — M^/M^) up to small corrections. To 
compute Af we first recall the general form of the electromagnetic vertex with Dirac and 
Pauli form factors [78j, namely, r^{k + q,k) = 7^Fi(g^) -|- ia^"qyF2{q^) /2M , where in our 
case Fi = and F2 = n. Using the conventions of Ref. [78], we introduce the polarization 
tensor 
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that 



With Af'"''" < 0.0010 at 95% CL [77], we find, e.g., with M = m. 
\k\ < 4.1 ■ 10"^ whereas if M = m^/lO that \k\ < 3.4 ■ 10"^ [H]. 

A variety of astrophysical constraints exist on the magnetic moment of the neutrino, and 
they can be adapted to our current case as well. They emerge, in particular, from the impact 
of the additional cooling mechanism such would render on stellar evolution and lifetimes 
and on supernovae [SDIEI]- An additional, albeit somewhat weaker, constraint comes from 
confronting element abundances with the predictions of big-bang nucleosynthesis, to yield, 
e.g., fiiy < 2.9-10~^°/iB [82]- These constraints can be significantly weakened by the candidate 
particle's mass [71]; the ability to produce particles of (9(10 keV) in mass and more in plasma 
at stellar temperatures is limited. In the case of big-bang nucleosynthesis, the constraints on 
the magnetic moment of a massive tau neutrino are also weakened, though values as large 
as /ij, ~ 10~^ fiB are nevertheless excluded [82] . 
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The constraints we have considered in this section can be weakened by other means as 
well. Since they arise from the effects of particle production, the most economical mechanism 
is compositeness; to include this, we need only include a form factor at each electromagnetic 
vertex. Thus we replace k — >■ — M|/M^)^, where is the compositeness scale, 

in our earlier formulae. Thus for M = mg/lO our earlier bound of |k| < 3.4 ■ 10^^ [IS] 
from Af relaxes to |k| < 1.5 if Mc = 2 GeV. In this scenario, however, the electrically 
charged constituents may well give rise to other observable effects. One possibility which 
avoids this would be to give a known charged particle, such as an electron, a small hidden 
sector interaction, so that it can help constitute dark matter, though its contribution to Af 
has already been taken into account. We proceed to consider the manner in which direct 
constraints can be set on fi. 



IV. GYROMAGNETIC FARADAY EFFECT 

A medium of free electric charges in an external magnetic field is circularly birefringent 
and gives rise to a Faraday effect [83] , as long as the light doe not propagate at right angles to 
the magnetic field. This effect has long been used in radio astronomy to study the properties 
of the interstellar medium [HI]. A Faraday effect can also arise in an electrically neutral 
medium in an external magnetic field, if the constituents carry magnetic moments and if 
they are aligned by that magnetic field to give the medium some net magnetization [201 EI] ■ 
We consider the latter possibility exclusively. 

To derive the gyromagnetic Faraday effect, we apply a magnetic induction in a mag- 
netizable medium with circularly polarized electromagnetic waves propagating parallel to 
it. The external field induces a magnetization A^tot, i-e., a net magnetic moment /volume, 
where A^tot = A4o + Ad and Aio results from Bq alone. The total magnetization of a 
medium at rest in the laboratory frame obeys the Larmor precession formula 

—J— = — ^A4tot X Btot , (6) 
at h 

so that g^M/^i noting = eh/2M with e > for a particle of mass M, is the gyromag- 
netic ratio of the magnetic-moment-carrying particle. We note /i = Sgfi^j, where S is the 
spin of the particle. The gyromagnetic Faraday effect was first derived for a ferromagnetic 
material [201 EI], for which use of the magnetic field H is appropriate. Since dark matter is 
only weakly self-interacting at most, our hypothesized dark matter should be treated as a 
paramagnetic material — so that we employ the magnetic induction B throughout, though 
the use of H is also commonplace [85]. Corrections to the Larmor formula result if the 
medium's particles move at a significant fraction of the speed of light, or if the particles 
possess a non-zero electric dipole moment [86]. We shall neglect the latter possibility and, 
moreover, shall consider dark-matter candidates for which relativistic effects are ultimately 
small corrections. 

To determine the relativistic corrections to Eq. ([6]), we first construct the covariant clas- 
sical equation of motion for a single spin in homogeneous electromagnetic fields. This is 
germane as we can and indeed do neglect the mutual interactions of the dark-matter parti- 
cles, so that the magnetization is given by the quantum-mechanical expectation value of the 
spin operator for a single particle times the number density; and the time evolution of the 
expectation value is itself described by that of the associated classical equation of motion. 
The latter, for a charged particle with a spin, is given by the Bargmann-Michel-Telegdi 
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equation [86]. We cannot use this result directly because no forces act on electrically neu- 
tral particles in homogeneous electromagnetic fields, so that no Thomas precession term is 
present [HZ|- Nevertheless, well-known treatments [83] can be readily adapted to this case. 
Requiring dW^ /dr = 0, where U = 07(1, /3) is the 4-velocity of the particle in the laboratory 
frame and S is its spin, namely S = (5°,«S), we find that 

dS- fF»% + ^{S,F^^U,)) , (7) 



where F"^ is the field-strength tensor in SI units and r is the proper time of the particle. 
Thus in the laboratory frame the magnetization evolves as 

dM-tot 



= ^ ( M,ot X {B,,, -(3x^)+ X Mtot) • (Stot - /3 X ^) ) , (8) 



where we emphasize 7 is the Lorentz factor, namely 7 = 1/a/1 — To proceed, we separate 
-Btot and Etot as -Btot = Bq + B and Etot = Eq + E, so that Aio results exclusively from the 
external electromagnetic fields Bq and Eq. Working to leading order in the small quantities 
Ad, B, and E, which arise in the presence of electromagnetic radiation, we have 

7^ = ^f>tox(So-/3x^) + Mx(Bo-/3x^)+Mox(S-/3x^) 
dt n \ c c c 

+7^/3 X Mo) ■ (Stot -(3x^) + {f3xM)-{Bo-(3x^)^y (9) 

We can consider the evolution of the dark matter magnetization in vacuum, i.e., in the 
absence of ordinary matter, or in matter. Since the largest external fields we can apply obey 
Bq':^ Eq/c in vacuum, we set Eq = henceforth. We note, however, that the atomic-scale 
separation of electric charges in matter permit the opposite limit, Bq ^ Eq/c, so that the 
analysis of the magnetization in that case can be altogether distinct. We set this possibility 
aside for later discussion and continue with the analysis in vacuum. We note that the ability 
to establish a vacuum relies on the presence of matter with conventional electromagnetic 
and strong couplings; dark matter is sufficiently weakly interacting that vacuum technology 
does not affect it. Hence we use "vacuum" to connote the absence of ordinary matter. In 
this case we apply Bq in the ^-direction and choose f3 to be parallel or antiparallel to x as 
well. As we have mentioned, the entry of dark matter into the magnetic field region acts as 
a spin filter device. The dark matter which does enter the apparatus can thus possess a net 
magnetization, so that AIq is in the i-direction^. As a result Eq. (|9| reduces to 

7^ = ^(>txBo + Atox(B-/3x^)) . (10) 

Choosing the wave vector k of the light in the ^-direction as well, we recall E = —cx x B 
and let B{x,t) = B±e±exp{ik±x — iurt), where e± = y ± iz. We define the polarization 



^ This follows irrespective of the sign of g. For g<0, however, the spins preferentially point in the ±x 
direction. 
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state with positive helicity, e+, to be right-handed, which differs from the convention used 
in optics. In steady state, we find A4. = A4±e± and finally that 



M± = ± B± = x±B± , (11) 

"fUJ ± UJb 

where we have chosen jS = —jSx and defined ujm = df^M-Mo/h and ujb = gfJ'MBo/h. Since 
k± = {uj/c)y/l + x±, "we have 



or that 



,^^0^1+^^ (12) 

where we note for conceivable light sources that uj ^ ujb,ujm- Physically, the magnetic 
field associated with the passing electromagnetic wave tugs on the spinning particle in a 
direction perpendicular to A^O) prompting it to emit radiation which interferes with the 
light traveling in the forward direction, generating the birefringence. The Faraday rotation 
angle (p is simply = (fc+ — fc_)//2, where / is the total distance travelled by the photon. 
The quantity um is signed, so that the sense of the rotation angle determines the sign of 
g. If (3 = 0, we recover the result of Ref. [12], whereas in the extreme relativistic limit, i.e., 
as 7 ^ oo, x± ^ c'-iid thus (p vanishes as well. The average value of k± is not altered to 
leading order in small quantities, namely, 

....4(...M^^(i + o(^)). (14) 

so that an appreciable Faraday rotation can accrue in the absence of an effect on the average 
group velocity. Thus far we have considered the Faraday rotation of linearly polarized light 
consequent to passage a distance / through a medium; practical considerations demand that 
we determine its properties under reflection as well. If we reverse the direction of the light. 



Eq. (11) is unaltered save for the sign of the term in (3. The last does change sign since 
E = —ck X B. Thus if we set /5 = 0, an initially right-handed circularly polarized wave, e.g., 
travels both forward and backward with wave number /c_|_, so that the rotation angle accrues 
coherently under momentum reversal. The additional Faraday rotation associated with the 



explicit /5-dependent term in Eq. (12), however, cancels under a round-trip transit. Such 
contrasting behavior is long familiar from the study of birefringence in chiral media |22], 
which break macroscopic parity invariance. Similar conclusions have been drawn from an 
analysis of parity- violating photon-external-field interactions as well [HE]- Thus the net 
rotation angle after a round-trip, or after many, of a total travel length / is 



l + Ohi^ ■ (15) 



^„=fli(l + 0(^, 

Z7C \ \7^co'"' 

Neglecting the O^u^"^) corrections and working in the /3 — > limit, this becomes simply 
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which agrees with the result of the non-relativistic treatment in Ref. [19]. For a system 
at rest in thermal equilibrium, the magnetization Ai^ is a simple function of the applied 
magnetic field. We recall that for a system of spin 1/2 particles, e.g., each with magnetic 
moment /i, the magnetization for a system with number density um at temperature T is 



A^o = ^MAitanh ^1^^ , (17) 

though it is of little practical relevance to the current circumstance, for the system we con- 
sider approaches thermal equilibrium extremely slowly. That is, a dilute gas in an external 
magnetic field polarizes through spontaneous emission, and the rate W for this process is 
given by that of a magnetic dipole transition [90]: 

4 ^ 2 



where uj = gfiMB/h. Thus even if fiM = /^b, the Bohr magneton, and i? = 25 T, we would 
have ly ~ 1 ■ 10~^ s~^, which is trivial compared to the average rate with which dark matter 
is expected to transit an experimental apparatus. Although dark matter may possess some 
primordial magnetization, it is likely so small [19] that it is important to realize other means 
of polarizing it. For the particular geometry we consider, as we have noted, the onset of 
the magnetic field region acts as a spin filter device. Although this method should yield 
some net magnetization for any non-zero spin S, we explicitly assume a spin 1/2 candidate 
in what follows. If the velocity vm of the incoming particles is aligned with the direction 
of the magnetic field, then particles with vm = Pc < fstop can only enter the magnetic field 
region if their magnetic moment is aligned with it, where fstop is such that 



1 

2^ 



Mv'=\fi\Bo. (19) 



We have set any external electric field to zero and have neglected corrections of 0{j3'^). The 
dark matter which does enter the apparatus can thus possess a net magnetization; namely, 

Mo = umIjV , (20) 

where V is the polarization of the spins. We define V = (A^+ — A^_)/(A^+ + iV_), where 
and are the number of spins pointing in and against the direction of Bq, respectively. If 



/i —fi then Aio — >■ A4o just as in Eq. (17). We study the value of P as a function /x-Bq, 
M, and astrophysical parameters in the next section. 

Before proceeding, we return to the notion of studying the Faraday rotation of dark matter 
passing through ordinary matter. In this regard, we wish to consider matter comprised of 
atoms with closed electron shells, so that there are no unpaired electrons present to engender 
a gyromagnetic Faraday effect. The exceptionally large electric fields associated with atoms 
and nuclei [91] ^, make it possible for \f3 x Eo\/c to exceed presently achievable external 
magnetic fields [92]. Such considerations yield significant limits on the neutron electric 
dipole moment [93], e.g., from neutron-noble gas scattering [HHIM]. Returning to Eq. (|9]), 



we choose x \\ k as in previous case, but now choose /3 _L i so that f3 x Eq can also be in 



In H-atom, e.g., a test charge a Bohr radius away from the proton sees E/c w 2 • 10^ T. 
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the X direction. Counting \ j3 x Eq\/c as a parameter of 0{1) and neglecting terms of 
and higher, we have 

dM g,uf (21) 



dt h \ \ c / V c 

and the steady-state solution 



= ±^^S± , (22) 



where ue = 5'/iA/|/3 x ^^qI/^c, which yields the rotation angle 

irrespective of whether round-trip paths are executed by the light. In this context, then, the 
effective magnetic field is relevant simply to the value of Aio- Here, too, we need to determine 
the polarization of the dark matter which penetrates the material. Suppose f3 = jSy and 
that it is possible to choose a material for which E^z ^ -Eqx, so that the effective magnetic 
field is in the x direction. The force on the dark-matter particle in entering the medium is 
V(^ • (/3 X Eq/c)); since we can expect the magnitude of Eq to depend on y, a longitudinal 
Stern-Gerlach effect is still possible in this geometry. A force in the z direction can engender 
the more familiar transverse Stern-Gerlach effect, but the increasing diameter of the laser 
beam as a result of scattering in its passage through the material may make it impossible 
to exploit this feature. Thus we tentatively conclude that it ought be possible to realize a 
meaningful Faraday rotation study in matter as well, with potential gains in sensitivity to 
a possible dark matter magnetic moment. We now return to the vacuum case, to describe 
how such an experiment can be realized and to estimate the limits on /z one could possibly 
obtain. 



V. A FARADAY ROTATION EXPERIMENT 

Although the Faraday rotation effect we discuss can be found through correlation studies 
of the polarization of the cosmic microwave background [TU], a terrestrial Faraday rotation 
experiment offers a number of advantages. Current bounds on primordial magnetic fields [39] 
make any primordial magnetization associated with dark matter small [19] and difficult to 
probe, particularly in experiments executed over terrestrial length scales. However, in this 
case, as we shall demonstrate, one can apply a very strong magnetic field of known strength, 
and polarize the dark matter to an appreciable degree. Moreover, as we have seen, the 
Faraday rotation associated with a magnetic moment accrues coherently under momentum 
reversal, so that the measurement can be made in a small cavity and yet have a long effective 
path length. Measurements of very small rotation angles are also possible; the recent PVLAS 
experiment, for example, was able to achieve sensitivity to rotation angles of O{10^^) [36] 
— this stands in constrast to a sensitivity of C(10~*^) anticipated with the future CMBpol 
satellite [95] . 

A schematic of the Faraday rotation experiment we propose is illustrated in Fig. [T} Its 
ingredients comprise a laser, a linear polarizer, an evacuated optical cavity through which 
the light makes multiple passes and to which a longitudinal, steady magnetic field has 
been applied, and an analyzer. The technical requirements of a sensitive Faraday rotation 
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FIG. 1: Schematic of the Faraday rotation experiment described in text. 



measurement may demand a more sophisticated setup; however, such details are not needed 
for our estimate of the hmit on /i for a given sensitivity to the rotation angle. We emphasize 
that since dark matter carries neither electric charge nor suffers strong interactions [9], it is 
unaffected by vacuum pumps and, indeed, can pass into cavities free from ordinary matter. 
The generic setup we propose, save for the nature of the magnetic field, is common to the 
PVLAS experiment [3S], which investigated the optical properties of the vacuum, as well. 
As that experiment was sensitive to extremely small changes in the photon polarization, we 
adopt it as a reference — we use certain of the parameters chosen in that experiment in 
order to estimate the achievable bounds on the magnetic moment as a function of the mass 
of the dark-matter candidate. A remark concerning the orientation of the apparatus in our 
schematic is in order. In our derivation of Eq. (16) we chose f3 = x, though our result is 
of more general validity. That is, if we return to Eq. (9^, setting Eq = and working to 
we see that the absence of the 0{l3) term in Eq. (16) results from the transformation 
properties of E under momentum reversal — regardless of the direction of /3. Consequently, 
we need not orient the direction of the magnetic field with respect to the dark matter "wind" 
in any particular way, as we indicate in Fig. [l| It will turn out, however, that the choice of 
the orientation of the apparatus with respect to the Earth's velocity can modify the efficacy 
of the longitudinal Stern-Gerlach device. 

In order to evaluate the sensitivity of the scheme we suggest to a dark matter magnetic 
moment, we need to make assumptions concerning its local velocity distribution and number 



density. Both of the latter quantities are relevant to the computation of Aio via Eq. (20) 



which gives rise to the rotation angle in Eq. (16). To do this we adopt the "canonical 
model" [nSl EZ] employed in the analysis of direct detection experiments [211 UBi ES] for cold 
dark matter. That is, we assume that the dark matter in our galaxy resides in a non-rotating 
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halo and that the velocity distribution function f{v) in that halo is that of an isothermal 
sphere. This assumption is motivated by simplicity. More realistic distributions can affect 
the expected event rates, as well as their temporal variation [HHl EHl HOP] . We note that 
V = Vm + where Vm is the velocity of dark matter relative to the Earth, which we 
introduced earlier, and v^ 'i^ the velocity of the Earth relative to the nonrotating halo of the 
galaxy. The value of v can range up to the galactic escape velocity, roughly 650km/s [96], 
though we do not impose this cut-off in what follows, because the normalization of the 
resulting distribution differs by less than 1% [97j. Thus the form of / is that of a Maxwell- 
Boltzmann distribution: 



The velocity vq is related to the root-mean-square velocity of the distribution; for a galaxy 
with a flat rotation curve it is argued to be equal to the radial velocity of the galactic 
disk |101j and thus in practice is taken to be vq ~ 220km/s [96|. The velocity i^s is 
determined by the sum of 



where is the velocity of the galactic disk in an inertial reference frame, Ug is the velocity of 
the Sun with respect to the galactic disk, and Uf. is the velocity of the Earth about the Sun. 
Assuming the Milky Way is axisymmetric, is fixed by the circular velocity at the Sun's 
radius from the galactic center, which is 220 ± 20 km/s [102^ I103j . In galactic coordinates 
(— f,f, z) |104j . in km/s, we thus employ [96]: Ur = (0,220,0), which follows if the Milky 
Way is axisymmetric, and Us = (9, 12, 7). The approximate speed of the Earth about the 
Sun is 30 km/s, whereas its approximate speed about its axis is 0.5 km/s. We ignore the 
effects of the Earth's rotation about its axis in our analysis, as it is no larger than the error 
in Us jl04l 1105] . Moreover, Ve is dominated by motion along the longitudinal coordinate 
I. Noting this and that the ecliptic is oriented at an angle of roughly 60° with respect to 
the galactic equatorial plane means we can approximate the piece of ue in the longitudinal 
direction fas 30 cos60° cos(27r((t - 152.5)/365.25)), to estimate ve ^ (232 + 15 cos(27r((t - 
152.5)/365.25)))i [96|. We retain this approximation for simplicity, though the presence of 
the other components of Ve, as well as a more realistic velocity distribution, are important 
to an assessment of the event rates, the size and phase of any temporal variations therein, 
and dark-matter exclusion limits to better than 0(100%) [HSl ESI 1106] . We note, too, that 
more precise determinations of the astronomical inputs also exist |104[ 1105] . 

Finally, to complete our description of the model, we choose a dark matter mass density 
of p = 0.3 GeV/cm [96]. The currently accepted range for p is 0.2-0.4GeV/cm'' [TOTlfTOS] 
for a smooth matter distribution with a spherical halo. Models of galaxy formation which 
relax the smoothness assumption can give rise to local dark matter densities both larger 
and smaller than this range |109] . where we refer to Ref. jllOj for a succinct summary 
of the possibilities. Some of the uncertainties in the model we outline can be correlated. 
For example, an elliptical halo, and concomitant triaxial velocity distribution, can yield 
somewhat higher local densities |11H 1112] . Unfortunately, direct information on the dark- 
matter mass density in our solar system is sparse, and observational bounds exceed the 
estimate we employ by orders of magnitude |113l 11141 11151 11161 11171 1118] . 

To realize a limit on /i, we assert in what follows that dark matter is comprised of a single 
type of spin 1/2 particle with fixed M and n, though we emphasize that the detection of 
a signal does not require that the particle have spin 1/2. Our numerical estimate is in two 




(24) 



Ve = Ur + Us + U, 



(25) 
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FIG. 2: The polarization fraction \ V\ after Eq. (27) plotted as a function of day for various fstop and 
vq. Here ve has been assumed parallel to Sq; this yields a geometry-independent lower bound to 
the polarization, as described in text. The solid lines have vq = 220 km/s and Vstop = 1,5,10,15,20, 
and 25 km/s, respectively, as one moves in the direction of increasing \V\. The dashed curves have 
^^stop = 5 km/s and vq = 200 km/s and vq = 240 km/s, below and above, respectively, the solid 
line with the same fstop- 



distinct parts. Operating in the canonical model, we first assess the polarization of the dark 
matter in the magnetic field region as function of M and /x. With this in hand, we can then 
determine the limit on with M which follows from a given limit on the magnitude of the 
Faraday rotation angle. 

We now proceed to estimate the polarization. We assume that the magnetic field is 
uniform in some direction x. If this is realized in a finite volume, with no magnetic field 
external to it, then a longitudinal Stern-Gerlach effect exists at each field boundary — i.e., 
the "wrong" spin state suffers a repulsive force at each surface. In what follows we assume 
a slab geometry and estimate the polarization resulting from crossing a single interface; we 
assume the magnetic field is perpendicular to the interface. In this case the polarization 
condition is on i>m ■ -Bo < "^stop, where fstop is fixed by Eq. ( 19 ). Thus the fraction of particles 
which enter the magnetic field region with 100% polarization is 



pol 



d?VM f{VM,VE)Q{v, 



stop 



\vm ■ Bq\) 



(26) 



so that its polarization in the magnetic field region is V = (sgn/i)/poi/(l — /poi), and its 
number density is hm = p(1 — fpoi)/M. The evaluation of Eq. (26) yields two pieces, which 
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are distinguished by whether \vm\ is larger or smaller than fstop- The term with \vm\ < "i^stop 
can be evaluated irrespective of the orientation of the magnetic field region, whereas the 



term with \vm\ > v^top depends on the value of Ve ■ Bq. With Eq. (24), we see that / 
is maximized for Vm in the neighborhood of —Ve, so that if fgtop ve, the term with 
I'^Ail > "^stop dominates /poi, and the relative orientation of Ve and Bq becomes important. 
If, moreover, we make ve -L Bq, then we can have both vm = —ve + S and vm ■ Bq = 6 ■ Bq 
with S small — we expect /poi to be maximized for this geometry. Thus to set a limit on /i 
irrespective of the orientation of i;^; and Bq, we choose Ve \\ Bq, as this geometry gives the 
smallest value of /poi for fixed /i and M. Ultimately the assessment of the limit on /x in a 
real experiment will depend on the geometry and orientation of the magnetic field, but our 
procedure should bound from above the limit on /i to be found from a given sensitivity to 
the Faraday rotation angle per unit length. In this, we implicitly assume that dark matter 
is described by a single constituent. In this special case we find 



Since erf{z) — >• 1 as z — *• cxd, /poi 1/2 as Vstop — > oo, as required. The form of Eq. (27) 
emerges from a partial cancellation of the contributions from the vm < ^^stop and vm > Vgtop 
regimes, so that we pause to consider whether the inclusion of an escape velocity in this 



particular case could modify our results. In this event, the normalization of Eq.(26) changes 
slightly, but negligibly [97], and the integrand accrues a factor of 0(fesc — \vm + ve\)- For 
^stop Ve, this additional factor does not restrict the region of integration unless vm ^ 600 



km/s; finally, we conclude that the continued neglect of Vesc is justified. Solving Eq. (19), we 
find Vstop ~ 4.51km/s a/ kBq [T] (rrie/M) , where Bq[T] is in tesla and nie is the electron mass. 
Employing a 7 T magnet, as used in the UCNA experiment |119] . though magnets of up to 
20 T are commerically available |120] . we find for M = rrie and k = 1, with fi = KfiM, that 



"Wstop ~ 6 km/s. The polarization fraction \V\ as a function of day, using /poi from Eq. (27), 



is shown for various f stop and vq in Fig. |2| Day-by-day variations in the polarization exist 
since /poi grows larger as \ve\ decreases. The time variation is more marked for the Ve \\ Bq 
geometry we have chosen, as the precise value of Ve impacts the range of velocities which 
enter the magnetic field region, and, hence, the polarization. For definitenss, we choose day 
335 to set limits on /j,. The experiment we consider can be sensitive to both annual and daily 
signal variations, in principle. This follows from the duration of the photon-dark-matter 
interrogation time, which, in turn, is set by the size and finesse of the cavity in which the 
experiment is realized. We note that the total travel length / of the laser light in the PVLAS 
experiment is / = 4.4-10^ cm [36], which corresponds to an interrogation time of (9(0.1 msec). 
This admits the possibility of using pulsed magnetic fields, which can be much stronger — 
as much as 89 T |121j . Irrespective of this, studies with the magnetic field on and then off 
are important to establishing a non-zero signal, if it is present. 

We now consider how a limit on the magnetic moment follows from our determined 
polarization and a given sensitivity to the Faraday rotation angle per unit length. For a 



spin-1/2 candidate, the rotation angle of Eq. (15) can be written as 



Vk' ^ (6.84 ■ 10-'\m')n^[cm-']l[cm]rK' Q)' , (28) 



with Ue = p(l — f pol) /"rrieC and V from /poi as per Eq. (27). We note, too, that if dark 



matter were a mixture of particle and antiparticles of the same mass that the rotation angle 
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FIG. 3: The limit on \n\ at 95% confidence interval, where ii = k^im and = eh/2M, as a 
function of M up to 1 MeV for various limits on (pQ /I and for different values of the magnetic field 
Bo. The solid lines correspond to i?o = 7 T, whereas the dashed lines correspond to Bo = 20 
T. The limit on (po/l in each case is lO^"*^^, 10^^'^, lO""*^^, and 10~"'^^ rad/m, respectively, at 95% 
confidence interval, as one sweeps from the top to the bottom of the figure. 



could cancel, at least in part. We can rewrite Eq. (28) in terms of a limit on by replacing 
Krrie/M with ^[^b\- The most stringent limits on \k\ emerge for M -C me- For fixed M, 
its numerical value rests on the ability to determine 0o/^- In the PVLAS experiment [36] , 
the error in (po/l is determined to be 0.5 ■ 10^^^ rad/m. Since / = 4.4 ■ 10^ m, 0o itself is 
determined to 2.2 • 10~^ rad. Assuming then that ipo/l can be determined to 1 ■ 10~^^ rad/m 
at 95% confidence interval, we find the limit on or as a function of M. In Fig. |3]we 
show the limit on |k| for candidate masses up to 1 MeV, whereas in Fig. |4]we show the limit 
on |/i| for masses from 1 to 100 MeV. The limits depend on the dark-matter polarization 
V as well. To illustrate the relative importance of V and the determination of (^o/^ to the 
limit on we show not only how the limits change if Bq is increased from 7 to 20 T but 
also, in Fig. [sj the value of |P| associated with each limiting value of \k\ in Fig. |3} The 
increase in Bq makes little difference at the lightest mass scales we consider, simply because 
the polarization at these scales is already near unity. Indeed, increasing the value of Bq 
simply increases the largest value of M for which we can reasonably constrain the value of 
\k\. Our assessment of the polarization as per Eq. (27) is that of a lower bound, yet the 
polarizations we find are large enough that our upper bounds on are no more than a 
factor of a few larger than what we would find after a realistic simulation of the geometry 
and orientation of the magnetic field region. It thus emerges that the limits to be set depend 
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FIG. 4: The limit on in units of fis at 95% confidence interval, as a function of M for masses 
ranging from 1 to 100 MeV for various limits on (po/l and for different values of the magnetic 
field -Bo- We employ the notation of Fig. [3] throughout, and note that dotted line corresponds to 
Bq = 89 T and a limit on (po/l of lO""*^^ rad/m at 95% confidence interval. 

overwhelmingly on the ability to determine 0o/^- Let us consider then the determination of 
this quantity and its consequences carefully. 

In Figs. ^ and ^ we show how the limits improve as the determination oi (po/l improves 
by orders of magnitude. This can be realized by either increasing / or by bettering the 
measurement of the rotation angle. We note, in particular, that better determinations of 
(po are possible |122j . and indeed that precision polarimetry at the shot-noise limit has been 
demonstrated |123l 1124] . In this limit the error in (pQ is set, crudely, by the number of 
photons counted, (50o — 1/(2 VJqT), where Jq is the number of photons per second and T is 
the measurement time. Assuming a 1 W laser in the optical regime, so that ^ 1 eV, we 
have 500 — 2 ■ 10~^° rad-s/ y'T[s] jl22l 1124] . Moreover, it has been demonstated that the use 
of squeezed light makes it possible to evade this quantum limit and realize yet more precise 
polarimetry |124] 1125] . All this suggests that the significant gains in the determination of 
(po/l we consider and more are indeed possible. Interestingly, as the determination of (po/l 
becomes more precise, increasing the value of Bq becomes more important to an improved 
limit on \k\. We note, e.g., that if M = 100 MeV and B = 89 T that the polarization 
fraction associated with the value of which follows from a limit of (po/l of 10~^^ rad/m 
at 95% confidence level is \V\ 0.09. 

We have determined the direct limits on or which would follow from a non- 
observation of Faraday rotation at a given sensitivity. As they stand the limits constrain 
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FIG. 5: The assessed polarization fraction \V\ associated with the Umiting value of \k\ for fixed 
M shown in Fig. |3j We assume, as in Fig. [2] that ve \\ Bq, to yield a geometry-independent lower 
bound to the polarization. The conditions which specify the various curves are as given in Fig. [Sj 

the possibility that the dark-matter particle is a composite built from electromagnetically 
charged constituents. Significantly more severe limits would be needed to challenge the 
possibility that a dark-magnetic magnetic moment exists by dint of quantum-loop effects. 
Thus we wish to consider the prospect of radical improvement in our limits. Improvements 



can come from bettering the determination of (f)o/l, increasing Bq, or, finally, noting Eq. (28), 
from increasing the number density Ue- The ability to limit the value of (po/l beyond that 
established by the P VLAS experiment [36] has already been demonstrated |123t 1124] , and the 
ultimate limit to the determination of (po/l has not yet been established |125] . Moreover, as 



discussed in Sec. IV , mounting the Faraday rotation experiment in matter can yield magnetic 
fields which are larger by orders of magnitude, and can possibly realize larger values of the 
dark-matter polarization V at fixed /i and M. Figures [3] and [5] give a sense of the outcome 
as a result of such improvements. Let us now turn to the last possibility. 

For the usually assumed Galactic halo density, ng is only some 600 cm~^, so that the 
number densities associated with warm dark matter are still very low indeed. We wish to 
consider whether the value of Ue can yield to experimental manipulation. In particular, we 
would like to increase the dark-matter number density in a particular region of momentum 
space. At first glance this would seem impossible |126j . because Liouville's theorem demands 
that the density of the phase-space fluid for a system governed by a Hamiltonian is a constant 
of motion. Yet the Maxwell-Boltzmann velocity distribution we assume, cum gravitational 
potential, is a solution of Liouville's equation, so that a modest increase of number density 
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can be realized by reducing the gravitational potential. Significant gains, however, are also 
possible, if inelastic processes operate, as in this case Liouville's theorem no longer limits 
the density. Such considerations are crucial to the construction of superthermal, ultra-cold- 
neutron sources |126] : however, adapting this technology to our current context does not 
appear to be practical. 

VI. SUMMARY 

A Faraday effect also exists for light transiting a dark medium of electrically neutral 
particles with non-zero magnetic moments in an external magnetic field [TH] . We have used 
this notion to describe a Faraday rotation experiment which can lead to the direct detection 
of dark matter were it to possess a magnetic moment fi. Alternatively, a null result can 
be used to limit the magnetic moment of a dark matter particle. The set-up involves an 
evacuated optical cavity to which an external magnetic field has been applied and through 
which the light makes multiple passes in the direction of the magnetic field. We assume the 
dark matter wind in the Earth rest frame sweeps through the apparatus. The size of the 
Faraday rotation angle is proportional to the magnetization of the dark matter, so that it is 
important to give the magnetic moments some net alignment. For particles with sufficiently 
low kinetic energy, the passage into the magnetic field region itself acts as a longitudinal 
Stern-Gerlach device; the highest energy state in the magnetic field can be barred from 
entering the apparatus, engendering a net polarization. We note that such a technique 
has been used to polarize ultra-cold neutrons with near 100% efficiency [23] . Employing the 
usual assumptions [96j concerning the mass density and velocity distribution of galactic dark 
matter employed in the analysis of existing direct detection experiments, we have estimated, 
for fixed astronomical input, the polarization of the dark matter in the magnetic field region 
as a function of its mass and magnetic moment, as well as of the applied magnetic field. 
Given a local dark- matter mass density, our limits on |/i| then follow from the strength of 
the applied magnetic field Bq and from the ability to measure (j)Q/l, the Faraday rotation 
angle accrued per unit length. 

We have studied a range of candidate masses compatible with dark matter as a warm 
thermal relic, namely, from 1 keV to 100 MeV. This window in candidate masses is not 
accessible via other direct detection techniques and thus is unique to our study. We find the 
strongest limits on emerge at the lightest mass scales we consider. In setting our limits 
we have employed the specifics of existing, related experiments as far as possible |231 ES], 
though we note that the possibility of more precise polarimetry has already been demon- 
strated [123^ I124j . The sensitivity of the limits we obtain are such that they constrain the 
possibility that dark matter is a composite with a magnetic moment, akin to a stable neutron 
without its strong interactions. Indeed this analogy has proven useful in adapting technology 
used to manipulate neutrons to our current case. The technical limits of the polarimetry 
measurements have not yet been established |124[ 1125] . and mounting the experiment in 
matter with concomitant gains in the applied magnetic field may prove feasible, so that we 
can ultimately expect better limits on than those found in this paper. 
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